Optical phonon scattering and theory of magneto-polarons 
in a quantum cascade laser in a strong magnetic field 
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We report a theoretical study ol the carrier relaxation in a quantum cascade laser (QCL) subjected 
to a strong magnetic field. Both the alloy (GalnAs) disorder effects and the Frohlich interaction are 
taken into account when the electron energy differences are tuned to the longitudinal optical (LO) 
phonon energy. In the weak electron-phonon coupling regime, a Fermi's golden rule computation 
of LO phonon scattering rates shows a very fast non-radiative relaxation channel for the alloy 
broadened Landau levels (LL's). In the strong electron-phonon coupling regime, we use a magneto- 
polaron formalism and compute the electron survival probabilities in the upper LL's with including 
increasing numbers of LO phonon modes for a large number of alloy disorder configurations. Our 
results predict a nonexponential decay of the upper level population once electrons are injected in 
this state. 

PACS numbers: 



I. INTRODUCTION 

A quantum cascade laser (QCL) is a semiconductor 
laser with specially designed cascade quantum structure 
to realize unipolar carrier transport and intersubband op- 
tical transition, and thus to achieve long-wavelength las- 
ing by overcoming the bandgap limit of the material 1]. 
Since the first observation of the population inversion in 
a QCL i, the mid-infrared, far-infrared and even THz 
QCL's have been realized though few of them operate 
at room temperature fS\. It has been observed that the 
performances of the QCL's deteriorate with temperature 
quickly when the emission wave length gets longer, a sig- 
nature of an increasingly detrimental thermal activation 
of non-radiative losses. Therefore it is important to un- 
derstand and control the non-radiative paths from the 
upper state of the lasing transition. A magnetic field 
B applied parallel to the growth axis has been used as 
an external parameter to monitor the characteristic of a 
QCL. In particular, it modulates the output power of the 
laser with an oscillation period proportional to 1 /B. This 
output is sensitive to both elastic and inelastic scatter- 
ings, which introduce non-radiative relaxation channels 
to carriers in the Landau level (LL) 0, H, H, 0, B Hi- 
Thus it is of interest to study theoretically the relax- 
ation mechanisms of the QCL in the presence of a quan- 
tizing magnetic field. Calculations of the LO phonon 
emission rates in these disordered systems have rarely 
been reported [13, [O, [13 • Yet, in the Gao.47Ino.53As 
active-region material the alloy effect nearly reaches its 
maximum value, with broadening of the LL up to several 
meVs. Therefore, it is worth devoting more efforts to the 
understanding whether the notion of phonon emission re- 
mains meaningful in LL quantized material. 

Assuming a LL structure at high magnetic field, the 
problem could be tackled in two radically different ways. 
In the weak electron-phonon coupling case, the electron 



states are initially computed by including alloy scatter- 
ing potentials and the phonon emission rates are then 
calculated using the Fermi's golden rule. This traditional 
perturbative computation would give an estimate of the 
exponential decay rate of the upper level electron popula- 
tion. On the contrary, LL electrons and LO phonons form 
magneto-polaron states when their coupling is strong 
enough. Similar effects appear recurrently in the opti- 
cal properties of quantum dots (QD's), e.g. [ll,[i3]- In 
this mixed mode description the electron-LO-phonon in- 
teraction is initially taken into account exactly. Hence, 
the very notion of phonon scattering/emission becomes 
irrelevant. We shall examine the magneto-polaron case 
where the alloy scattering broadens the magneto-polaron 
states and compute the time-dependent survival proba- 
bilities in the upper LL's for an electron initially in one 
of the alloy broadened upper LL's. This electron-phonon 
strong coupling theory leads to a distinguishably different 
prediction for the carrier relaxation process as composed 
to the ordinary phonon emission picture, and predicts a 
nonexponential population decay from the upper levels. 



II. LO PHONON EMISSION OF THE ALLOY 
BROADENED LL'S 



For mid-infrared QCL's, it is a reasonable assumption 
that the alloy scattering destroys neither the subband 
structure at zero field nor the LL structure at high mag- 
netic fields. So we consider only a few subbands and their 
LL's. Leuliet et al demonstrated that the alloy scatter- 
ing effect in the barrier is negligible as compared to that 
in the well ll]. We correspondingly model the active 
region of a QCL as a single quantum well (QW) clad 
between infinite potential barriers along the z direction. 
A uniform magnetic field applied parallel to the quan- 
tum confinement direction is considered in Landau gauge 



2 



(with the vector potential A = Bxy). Moreover, we ne- 
glect the spin Zeeman effect. Thus, in the effective mass 
model the unperturbed degenerate LL eigenstates of the 
Zth subband read 



{r\Ei,n,ky) = -^=e'''yyxi{z)ipn{x + X^ky), 



£l,n — El + 



(1) 



where xi the wavefunction of the Ith subband, 
(/?„ the nth Hermite function, A — \J h/eB the mag- 
netic length and Wc = eB /m* the cyclotron frequency. 
We take m* — O.OSmo for electron effective mass and 
Lz = 12.33nr7i for quantum well width, which gives the 
zero-field subband separation £^2 — Ei = 147.5meV. The 
LL degeneracy (per spin) is D ~ L^Lj,/(27rA^). 

For a numerical calculation of the al- 
loy scattering effect in the Gai-xIuxAs well, 
we partition the system (a large box of 
no = X Ly X ^ 99.79nTO x 99.79nm x 12.33nm) 
into tiny unit cells (wq = 2.93Ax2.93Ax5.87A) 
In each cell the scattering potential is a random vari- 
able which equals to xAy with probability 1— x and 
to — (1— x)Ay with probability x, where x= 0.53 for 
Gao.47Ino.53As and AV — O.GeV On the scale of 

LL's or z dependent wavefunction extensions, the alloy 
fluctuations act like delta scatterers and we assume 
there is no correlation between the x values of different 
cells. Thus the alloy potential is given by [l7[ 



Valloy = AVUJO ( J2 ^'^(^- ^Ga) - ^ (1 - x)5(f - 



(2) 

With a given configuration of the alloy disorder, we cal- 
culate the alloy broadened LL's by diagonalizing the alloy 
potential in the basis of the degenerate LL's (for details, 
see Ref. [iBl)- With our sample box this means [D] = 59 
states per LL at 24.58T. The diagonalization has been 
done with N = 100 samples with randomly distributed 
alloy atoms. 

The numerical results (Figs. 1 and 2) demonstrate the 
validity of the LL structure and clearly show the localiza- 
tion of the electrons at the broadened LL's tails as well as 
the extended behavior at the maximum of the density of 
states (DOS). Fig. 1 shows the in-plane electron density 
distributions P{p) for four states of one particular alloy 
configuration, as defined by 



Pi,nAP) = / Xl{z)dz\^|Jn^^{p}\ 



(3) 



where (r|i5;, rt, i/) — Xi[z)''l'n,v{p) is the vih disordered 
state related to \Ei,n) (a short notation for all unper- 
turbed or disordered states related to the Zth subband 
and nth LL's). For each of the n, z/) states of a par- 
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FIG. 1: (Color online) Typical in-plane electron density dis- 
tributions (in A~^) for four states in the set of independently 
broadened LL's at 24.58T: (a) DOS tail of 1^2, 0); (b) DOS 
max of 1^2,0); (c) DOS tail of j£i,2); (d) DOS max of 
|Si,2). 
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ticular run, we have computed the x- and y-uncertainties, 
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FIG. 2: (Color online) Position uncertainties of the indepen- 
dently alloy broadened LL's at 24.58T: (a) 5^ of j-E'2,0); (b) 
5y of \E2,0)\ (c) 5, of \Ei,2) ; (d) 5y of \Ei,2). 
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Fig. 2 shows 6x and 6y for all the independently alloy 
disordered states related to |£^2,0) or |i?i,2) at 24.58T 
{[D] — 59; N = 100 disorder configurations). For an un- 
broadened LL |_E;,n, with the electron entirely de- 
localized along the y axis over a segment L there is 

Ay/ n + |. For L = 10^ nm and 

B = 24.58r, this means 5y — 28.8nm, S^"^ — 3.66nm and 
(2) 

6x — 8.18nm. Fig. 2 shows that the tail states are con- 
siderably more localized than the central states. Indeed, 
for the latter Sy is close to the expected value for a plane 

wave. We also note that ^i"-* and are considerably 
larger than 3.66nm and 8.18nm, respectively, and in fact 
close to Sy. This proves that the extended states of the 
alloy broadened LL's are very different from the unper- 
turbed \Ei,n,ky') states but display similar extensions 
along the x and y directions (the small discrepancies are 
due to the different boundary conditions along the x and 
y directions). 

We consider the electron-phonon interaction here only 
due to the GaAs-like LO phonons, because in the 
Gao.47Ino.53As alloy the Frohlich interaction is dominant 
and the potential of GaAs-like modes dominates over that 
of the InAs-like modes [l^] . We take the zone-center LO 
mode frequency as Huilo = SSJmeV and the Frohlich 
interaction as 



(6) 



where bt is the creation operator of an LO phonon 
with a wavevector of q and the Frohlich factor is 
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Ml 



2eoVa, 



£s^) with £00 = 11.6 and Eg 



13.3 



as the high-frequency and static relative dielectric con- 
stants of GaAs material, respectively, and ¥„■ the crystal 
volume. 

In the absence of the alloy-disorder, the un- 
perturbed factorized states |i?2, "2 = 0, fc^) (g) |0lo) 
and Til = p 7^ 0, fc,y) (g) |1lo) cross at the field 

Bp = m*{E2 — El — huj Lo) / ip^f^) for any values of ky 
and k'y. Let us consider the p — 2 resonance, for which 
the turn off of QCL's lasing has been numerically cal- 
culated under the assumption of a dominant inhomo- 
geneous broadening [llj. We can calculate the scatter- 
ing rate due to the electron-phonon interaction from one 
of the alloy broadened \E2,0) LL states towards any of 
the I -El, 2) broadened states assuming the Fermi's golden 
rule holds. The scattering frequencies are averaged over 
different alloy configurations. As shown in Fig. 3, at 
resonance {B2 = 24.58T) the scattering rate of this irre- 
versible departure from the initial state reaches a maxi- 
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FIG. 3: LO phonon emission rate from alloy broadened LL's 
\E2,0) to alloy broadened LL's \Ei,2) at B2 = 24.58T, as a 
function of the energy of the initial state. 



mum value when the initial state coincides with the cen- 
ter of the broadened LL's; this energy also corresponds to 
the largest DOS of the final states. Note that in spite of 
the fact that we deal with transitions between different 
LL's associated with different subbands, the scattering 
rate is quite large: we find in Fig. 3 a maximum scatter- 
ing rate of 8.66ps~^. This relatively large value results 
from the increased DOS of final states due to the Lan- 
dau quantization even in the presence of disorder. It also 
points out a difficulty in applying blindly the Fermi's 
golden rule since the energy uncertainty h/r « 5.70meV 
due to LO phonon emission is comparable to the LL 
width due to alloy broadening (with the full width at 
I/e DOS max as A « 5.69meV, see Ref. and also 
Sec. III). 



III. ALLOY BROADENED 
MAGNETO-POLARON STATES AND THE 
SURVIVAL PROBABILITY IN THE UPPER LL'S 

A. Unperturbed magneto-polaron states 

As shown in the Fig. 4, the Frohlich inter- 
action replaces the crossing of |i?2, 0) (8) |0lo) and 
\Ei,p) (81 |1lo) by an anticrossing and leads to the for- 
mation of mixed magneto-polaron states. These are ob- 
tained by diagonalizing the Frohlich interaction in the 
factorized states manifold [l^, and read: 

\±p,ky) = a±p\E2,0,ky) <^\0) 

+ ^f3±p,q\Ei,p,ky - qy) 



P±p.q = -i ^""^^ I*{q), 



S / ^'^ 

£±p = £2,0 + fy^±p , u;±p = ^ ± y ^ + ^1, (7) 
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FIG. 4: Illustrations of magneto-polaron states in ideal 
systems. (a) A schematic picture of the formation of 
magneto-polaron states. Insert: energy anticrossing around 
B2 = 24.58r. (b) Polaron splitting energies as a function of 
Bp for p = 1 to 10 (the line is only a guide for the eyes) . 



where ± refers to the upper {+) and lower (— ) po- 
laron states in Figure 4(a), the detuning energy is 
hSp = Si^p + hujLO — £2,0, and the half polaron splitting 
energy is given by 



hnp^ jY,9^\Ip{q)\^/q\ 



(8) 



with the electron-phonon interaction coupling matrix el- 
ement as 

Ip{q) = {E2,0,ky\e'^-''\Ei,p,ky-qy) 
{\{iq^ - qy)Y 



-exp 



( 



iX qxiky — 



(9) 



Actually, after the diagonalization, we find several un- 
coupled LL's one-phonon replica states in addi- 
tion to the magneto-polaron states. These uncoupled 
states cannot be populated by electrons entering the ac- 
tive region directly from the zero-phonon states of the 
injection region. For simplification, in this section we 
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concentrate on the magneto-polaron states only (the role 
of uncoupled states is discussed in Sec. IV). The electron- 
phonon scattering obeys the momentum conservation 
law. Since each of the unperturbed LL's has a definite 
ky, the Frohlich interaction conserves the ky number and, 
consequently, each magneto-polaron state has a definite 
ky value in the ideal system. For the case of the exact 
resonance between |i?2,0)(g) \0lo) and |1lo) 
the upper and lower polaron branches |±p, fcj,) have 
equal weights of phonon components while for a finite 
energy detuning these weights are unequal. This is sim- 
ilar to the case of strong coupling between atoms and 
photons in cavity QED except that each polaron branch 
has a degeneracy of the LL's The polaron splitting 
energy has a Vb dependence and is about A.OmeV at 
B2 = 24.58T [ig|. 



B. Disordered polaron DOS: 
SCBA and numerical results 

The alloy disorder breaks the ky-ie\a.ted translational 
invariance and thus either broadens the magneto-polaron 
states (weak disorder) or even destroys the polaron 
picture (strong disorder). As for quantitative esti- 
mate, by using the self-consistent Born approximation 
(SCBA) [2l|, we solve a pair of coupled equations of the 
self-energy of polarons given by 



S+p(£) 



S-p(£) 



T/2 

' +P,+P 



^ +P--P 



e — e+p ■ 



S+p(e) 



S-p(e)^ 



(10) 



V? 

+p,—p 



e — e+p — Y,+p(e) e — e-p — S_p(e) 



(11) 



where the alloy squared scattering matrix element is av- 
eraged over the alloy fluctuations viz. 



alloy 



(12) 



For infinitely large (2-dimensional) samples we have 
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where V^ = x( 1 — x) ( A F ) and each of the alloy squared 
scattering matrix elements involves two alloy scatter- 
ing events and consequently contains three terms: zero- 
phonon term when both events occur among zero-phonon 
states; one-phonon term when one event occurs among 
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FIG. 5: (Color online) Alloy disorder broadening of magneto- 
polarons: (a) SCBA calculation of the polaron DOS for p — 1 
to 4 at the resonant fields Bp; (b) SCBA calculation of the 
polaron DOS (in unit of m*/{TTh^)) around B2 = 24.58r; (c) 
and (d) are numerical results at 23. SOT and 24.58T, respec- 
tively, to compare with the SCBA calculations. 
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zero-phonon states and the other among one-phonon 
states; two-phonon term when both events occur among 
one-phonon states. The one-phonon term has a common 
factor as 



9 



E 



2P-ip! ^ 



■{\q„) exp{~-\^q)), 



and the two-phonon term's factor reads 
3.9^ 



(15) 
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, fe \ 2 
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1 



'exp{- -X^Aqj) 



(16) 



Once the magneto-polaron self energy is obtained by 
iteration, the polaron DOS near Bp is given by 



D 

pie) — Im 

TT \e 



1 



1 



e+p — S+p(£) £ - 



S-p(£) 
(17) 

As shown in Fig. 5(a), the alloy broadening effect is 
enhanced with increasing magnetic field. For larger p the 
alloy broadening may be smaller than the polaron split- 
ting. However, when p is large enough Bp is very close to 
Bp+i so that the two-LL's formalism [Eqs. (10) and (11)] 
is no longer valid. To set a criterion, we find that a finite 

polaron gap will open when 



2Vi 



< h^„. We 



can see from Fig. 5(a) that this criterion is not fulfilled 
for p = 1 and only marginally for p = 2, but holds for 
larger p values with recognizable double maxima in the 
DOS. In Fig. 5(b-d) we show more results for the SCBA 
polaron DOS of p = 2. 

We have also performed a numerical diagonalization of 
the alloy Hamiltonian for a given realization of disorder 
within a basis of polaron |±, ky'^ states. Once the polaron 
levels are calculated the density of states is calculated and 
an average over N=100 samples is performed. As seen 
in Figs. 5(c) and 5(d), the numerical calculations agree 
quite well with the SCBA results though the broadening 
effects appear larger in the numerical outputs. We can 
see in both SCBA and numerical calculations that out 
of resonance (Fig. 5(c)) the DOS displays two peaks of 
uneven heights. 

The DOS peaks of the two broadened LL's appear very 
asymmetrical at large detunings. The polaron level that 
resembles the I-E2, 0) <S) \0lo) LL (far from resonance) ac- 
quires a width and a shape that goes smoothly to those 
of the pure electron |i?2,0) LL [1^. On the other hand, 
at large detuning the SCBA broadening becomes insuffi- 
cient to fully account for the alloy broadening. This is be- 
cause a large number of uncoupled states become nearly 
degenerate with the polaron levels that resemble the one- 
LO-phonon replica states \Ei,p') (g) |1lo) and these un- 
coupled states are not accounted for neither in the SCBA 
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FIG. 6: Survival probability Ps for B2 = 24.58r with the 
initial state as the central level of the broadened |i?2, O) LL's 
while neglecting all the uncoupled states. 



equations nor in the numerical results in Fig. 5 (see Sec. 



IV). 



C. Time-dependent survival probability 
in the upper LL's 

The irreversible emission of LO phonon is the most 
efficient energy relaxation mechanism in QCL at zero 
magnetic field provided it is energy allowed. In the pres- 
ence of a strong field, in an ideal sample, this irreversible 
emission is impossible because of the formation of inter- 
subband magneto-polarons. Alloy scattering blurs the 
magneto-polaron states and it is interesting to ascertain 
the nature of the phonon emission in QCL subjected to 
a strong quantizing magnetic field. To this end, in the 
case of the p = 2 resonance, we compute for a given re- 
alization of the alloy disorder the survival probability in 
the ensemble of independently broadened LL states |z/) 
related to |£'2,0) once the carrier has been injected in 
one particular state |i^o) of this ensemble: 



-j—j \vo) ® \0lo) 

(18) 

where H ^ He + Vaiioy + Hph + He-ph- In the following, 
we take [I'o) as the central level of the broadened ensem- 
ble. Once this is achieved we average over the N = 100 
realizations of the alloy disorder. Several magnetic fields 
have been considered. 

When there is a large energy detuning between the 
zero-LO-phonon replica of |i?2,0) LL's and the one-LO- 
phonon replica of = 2^ LL's, the survival probabil- 

ity very quickly oscillates between 1 (t = 0) and Pmm to 
stabilize to Pqo = (1 + Pmin)/2 in a characteristic time r 
of a few picoseconds: if |i/o) is the central state of the 
broadened LL's and for B = 20T, we find r « 2.8ps and 
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Poo « 0.83. Right at resonance B2 = 24.58T, see Fig. 6, 
the averaged survival probability displays a large number 
of fast oscillations (pseudo period: 1.03ps) with a decay- 
ing amplitude (t « A2ps) that brings P, to the limiting 
value Poo = 1/2. We note that the pseudo period is close 
to 7r/il2 = 1.04ps. The limiting value 1/2 recalls the fact 
that the Frohlich interaction being so efhcient that the 
carrier very quickly exits from the initial state. However, 
once in the one-LO-phonon replica of the p = 2 LL, there 
is no sink mechanism the electron can use to escape from 
this LL. Thus, another Frohlich interaction brings it back 
into the |i?2,0) LL. There is no reason a priori for this 
state to be the same as the initial state. Subsequently, 
the electron leaves the initial LL and oscillates back and 
forth between the zero-LO-phonon replica of |i?2, O) LL's 
and the one-LO-phonon replica of |£'i,p = 2) LL's. This 
oscillatory cycle between the two broadened LL's is how- 
ever irreversible to the extent that the initial state has 
very little chance to be recovered after 2j jumps. This 
is in striking contrast with the coherent polaron oscilla- 
tions that would result if there were no alloy scattering 
acting to blur the polaron states, thereby leading to an 
oscillatory cycle between only two polaron states. 

The consequence on our understanding of the exper- 
imentally observed QCL oscillatory output versus B is 
significant. Unless one can find a plausible and fast es- 
cape mechanism from the one-LO-phonon replica, there 
is no reason to invoke an irreversible escape from the up- 
per state of the lasing transition in a QCL. The long lived 
oscillations displayed in Fig. 6 will actually be limited 
by the finite lifetime of the phonons due to anharmonic- 
ity effects. This lifetime is about lOps at low tempera- 
ture and a few picoseconds at room temperature in bulk 
GaAs. Similar values were found in InAs self organized 
QD's [25|. Beyond this lifetime, on average, the oscilla- 
tions should stop and the electron should have relaxed to 
one of the disorder broadened LL states | i?i , p) with no 
phonon. Note that such a relaxation path would be in 
complete contrast with bulk and QW materials at zero 
magnetic field but much in agreement with the energy 
relaxation scenario established in QD's [l3 |. 



IV. EFFECTS OF THE UNCOUPLED 
ONE-PHONON STATES IN THE STRONG 
COUPLING REGIME 



The approximation of neglecting the uncoupled one- 
phonon states turns out unsatisfactory because of the 
very efficient alloy scattering between magneto-polarons 
and uncoupled one-phonon states. As a result, the Rabi 
oscillations due to polaron effects will be severely damped 
since the alloy scattering between polaron states and the 
large number of uncoupled one-phonon states will give 
fast additional leakage channels for the LL electrons. The 
uncoupled states related to the | ±p, /cj,) polaron are ap- 



proximately 

^uncoupled) 

^uncoupled 



\Ei,p, k. 



y 

hhJLO 



qy) ® \ ^q),ky 7^ ky, 



(19) 



There are Nuncoupied — i[D] ~ phonon such uncoupled 
states for each pair of \±p,ky') (where Nphonon > 480 
discretized LO phonon modes should be included to re- 
produce accurately the polaron gap, and an accurate 
numerical diagonalization would give the number in to- 
tal as Nuncoupied = [D]{Nphonon - 1))- Notc also that 
these states are energetically placed between the two po- 
laron levels I ± p, fcy): at their mid-distance at resonance 
B = Bp and tend towards the fast increasing (with field) 
polaron branch at high detuning | B — Bp | . 

The effect of alloy disorder is twofold: (i) broaden- 
ing the uncoupled level and (ii) admixing the uncou- 
pled with the polaron states. However, a full diag- 
onalization, including all polaron and uncoupled one- 
phonon states with a large number of phonon modes, 
is too heavy numerically (this would mean diagonaliz- 
ing for each run of the N = 100 a matrix of the size as 
Dim(i7) = {[D]{NpHonon + 1))' > (2.8 x 104)2). In the 
following, we present estimates illustrating the impor- 
tance of effects (i) and (ii). 

An estimate of the effect is obtained by remarking that 
at the lowest order in Vaiioy it produces a shift of the 
polaron states which reads 



S±p 



{ ^uncoupled \{ ± p\VaUoy\unCOUpled)\^) 



■Hoy 



^uncoupled 



(20) 

It is difficult to give a simplified analytical expression of 
s±p hke Eqs. (13-14) due to the large value of Nphonon re- 
quired. With the results of the numerical diagonalization 
as shown in Fig. 5(d), this shift is roughly OAmeV at res- 
onance B = B2, and the effect (ii) becomes very impor- 
tant at increasing detuning | hS2 1 to the polaron branch 
whose e polaron tcnds towards e uncoupled- For increasing 
I hS2 1 it becomes necessary to take into account higher 
orders' effects of Vaiioy, once the uncoupled one-phonon 
states enter the DOS of the closer broadened polaron 
branch. Thus it turns out to be extremely inefficient to 
handle the problem by making corrections to the SCBA 
Eqs. (10-11). 

Another way to evidence the importance of the cou- 
pling between the polaron and the unadmixed states con- 
sists in looking at the irreversible escape from one of the 
polaron states to the unadmixed states broadened by al- 
loy scattering. An estimate using Fermi's golden rule 
gives at resonance B = Bp 



,P,ky~ qy\Valloy\Ei,p, ky)\ 

xp{e = e±p) , (21) 



2ttt±p 
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FIG. 7; (Color online) Electron-LO-phonon coupling element 
related to Frohlich factor versus phonon wave vector. 



where the DOS for one of the unadmixed states is written 



Pis) 



^uncoupled 



(22) 



Putting numbers we find t±2 ~ O.Sps, i.e. comparable to 
the polaron period. Thus, it is invahd to neglect the un- 
coupled one phonon states and effect (ii) is strong enough 
to generate states without well-defined polaronic charac- 
ter. Note that this result is peculiar to the alloy scat- 
tering we have used for static scatterers. In GaAs based 
QCL's the alloy scattering is negligible and the interface 
defects are definitely weaker scatterers [13, El- Hence, 
for those materials the polaron levels should be signifi- 
cantly more long-lived. 

Therefore, we have included the uncoupled states 
in the numerical computation of the electron time- 
dependent survival probability in the upper LL's |i?2,0) 
(Eq. (18)). We take a certain number of phonon 
modes and ensure that the half polaron splitting energy 
fi^2 > l.SmeV since a small portion of the phonons 
give dominant contributions to the electron-phonon 
coupling (Fig. 7). As an approximation, we choose 
a smaller quantization size in the quantum well free 
dimension (viz. aL^ x aLy x L^) for the phonon modes 
than for the LL's electrons (which is Lx x Ly x L^). 
In the z-dimension, a specified wavevector is chosen 
for the phonon modes (typically 2-k/Lz). Numerical 
computations of the upper LL's survival probabilities 
with initially putting an electron in the center of the LL's 
I £^2,0) have been done with a = 1/5 (iVp^o„o„ = 20, 
Dim(i/) w (1.2 X 103)2), a = 1/4 {Nphonon = 36, 
Dim(i/) w (2.2 X 103)2), a = 1/3 {Nphonon = 56, 
Dim(H) « (3.4 X 103)2) and a = 1/2 {Nphonon = 88, 
mm{H) « (5.3 X 103)2) (Fig. 8(a)) for B2 = 24.58T 
(for a = 1 more than 384 modes are required). We see 
that at resonance (Fig. 8(a)) the survival probability 
in the |i?2,0) LL's decays very rapidly with time: in 
about 0.6ps it has dropped to nearly zero. The decay is 
almost insensitive to the number of uncoupled phonon 
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FIG. 8: (Color online) Taking increasing numbers of the un- 
coupled one-phonon states into account, computations of the 
electron survival probabilities with initial state in center of 
the alloy broadened LL's |i?2,0): (a) remaining in j_B2,0) at 
B2 = 24.58r; (b) in polaron states | ± ) at B2; (c) in |-E2,0) 
at 20.00T; (d) in |J52,0) at 23. SOT (see text). 
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FIG. 9: (Color online) Computations with 36 phonon modes 
of the long time limit of the survival probability Pi,(oo) 
in the |i52,0) LL's: (a) Initial states in center of LL's 
|-E2,0) ® |0lo) for each sample; (b) Average over all initial 
LL's |£2, 0) Cgi |0lo) for each sample. Also the ensemble aver- 
ages are taken over the N=100 samples for each case (jointed 
by lines, which are guides for the eyes). 

modes and in fact coincides with the one shown in Fig. 
6. However, instead of recovering to 1/2 at long time, 
Ps remains very small and displays faint remnants of the 
polaron oscillations. This is the irreversible oscillation 
between the polaron branches and the Nuncoupied uncou- 
pled phonon modes that explains the lack of recovering 
after the first oscillation. 

In Fig. 8(b) we show the survival probability of re- 
maining in the polaron states with the same initial state 
as that of Fig. 8(a). It should be noticed that at the first 
half picosecond the polaron effect is still strong enough 
to hold its pseudo oscillation against the damping due to 
the alloy scattering effect. The same calculation of the 
probabilities in |-E2,0) for lower magnetic fields near B2 
are shown in Figs. 8(c) and 8(d). The Figs, reveal that, 
out of the resonant field and in the presence of both alloy 
and LO phonon scatterings, the probability of relaxing to 
the lower LL's remains small at long time. 

The experimental implication of our findings on the 
resonant non-radiative relaxation in the QCL's subjected 
to a quantizing magnetic field can be illustrated in the 
following way. For one realization of disorder and a given 
magnetic field B, we compute the long time limit of the 



survival probability Ps(oo) in the |-B2,0) Landau level 
in two cases. In case (a), the initial state is in the cen- 
ter of the alloy broadened |£^2,0) ® |0to) LL's. In case 
(b), the quantity Ps(oo) is computed for every state of 
l-E'2, 0) (g) |0lo) and its average < Ps(oo) > over all these 
initial states is computed. Then another realization of 
the disorder is created and both Ps(oo) and < Ps{oo) > 
are calculated again. Fig. 9 shows the curves Ps(oo) and 
< ^5(00) > versus B. If alloy scattering and electron- 
LO-phonon interactions are the dominant non-radiative 
losses, the curves shown in Figs. 9(a) and 9(b) can be 
compared to the output power of the QCL. It appears 
that the calculated widths are comparable to those seen 
experimentally [Tol. ITlj. 

V. CONCLUSION 

We have presented a theoretical analysis of the 
electron-LO-phonon interaction in QCL's structures in 
the presence of a strong magnetic field and of static short 
ranged scatterers. Our objective was to ascertain the ac- 
curacy of the weak coupling regime between electrons and 
LO phonons. Our results show that the very notion of a 
Fermi's golden rule is highly questionable in these struc- 
tures because of the LL singular density of states. Neither 
the static scatterers are weak enough nor the electrons 
and LO phonons form stable magneto-polarons. As a 
result the survival probability in the upper state of the 
lasing transition never decays exponentially to zero but 
displays a number of damped oscillations before stabiliz- 
ing to 1/2, thereby evidencing that the magneto-polaron 
levels never completely empty. Or, the static scatterers 
are efficient enough to wash out the polaron oscillations 
because it couples these polaron levels more efficiently 
to the huge reservoir of uncoupled one-LO-phonon states 
than to the polaron states. The survival probability in 
the upper state of the lasing transition decay to zero but 
not at all in an exponential fashion as would result from 
the Fermi's golden rule. Instead, it first decays like in the 
damped polaron case and once has reached its first min- 
imum practically never recovers. The complicated time 
evolution of the survival probability evidences the need of 
a more microscopic description to understand the physics 
of the non-radiative mechanisms in QCL's. It also warns 
against the estimated efficiency of the static scatterers 
or phonon scattering when it is based on oversimplified 
descriptions of the disorder on the QCL quantum states. 
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